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ELECTRON EMISSION PRODUCED BY PHOTOINTERACTIONS 

IN A SLAB TARGET 

by Byron E. Thinger* and James A. Dayton, Jr. 

Lewis Research Center 

SUMMARY 

A method of estimating the production and penetration of energetic electrons in a 
target as a consequence of photon interactions is described. The calculations incorpo- 
rate the experimentally determined photoelectric, Compton, and pair production inter- 
action cross sections and utilize empirical electron energy range relations. The predic- 
tion of energetic electron emission is based on two quantities which are derived in this 
report, the probability of electron escape and the average path length of escaping elec- 
trons. 

In this analysis the fission gamma energy spectrum is divided into multiple energy 
groups, each of which is defined in terms of a monoenergetic gamma flux. With the 
mean energy of each gamma group as a reference, the number and average energy of 
primary electrons produced within the target are determined. The transport and loss of 
energy of those electrons which are able to travel through the target to the surface are 
computed to give the power per unit area of the emerging electrons and the correspond- 
ing escape current density. These latter two quantities are of interest in an application 
to thermionic conversion within a nuclear reactor system. 

The results are presented for a numerical computation of the current density and 
energy flux of the electrons emitted from a 0. 33-centimeter-thick tungsten target irra- 
diated in the core of the NASA Lewis Research Center Plum Brook nuclear test reactor. 

INTRODUCTION 

The ionizing effects of gamma radiation have been studied as a means of reducing 
the plasma impedance of in-core nuclear thermionic convertors by Butler (ref. 1) and 

* Assistant Professor of Engineering, California State University, San Francisco, 
California; Summer Faculty Fellow at the Lewis Research Center in 1969. 



Forman (ref. 2). Gas -filled electronic tubes appear to conduct more current when they 
are operated in or near a nuclear reactor. Presumably the conductivity of the filling 
gas is increased because of ionization of the neutral gas atoms by energetic photons and 
by electrons generated by the interaction of gamma rays with the tube surfaces. For- 
man's initial reactor experiments indicated some promise for thermionic converter per- 
formance improvement. However, Butler’s analysis predicts that the ionization con- 
tributed by the gamma radiation field is small, and no enhancement of operation could 
be expected. 

The primary objective of this work is to calculate the current density and energy 
spectra of electrons which are emitted from the surface of a metal target placed in the 
gamma radiation field in the core of a nuclear reactor. Since experiments within nu- 
clear reactors are quite costly and require extensive hazards analysis and lead time, 
other methods of studying the effects of gamma radiation on the performance of thermi- 
onic converters are of interest. For instance, one method might be to simulate the yield 
of primary electrons from a target by using the beam from a linear accelerator, for 
example, the 3-Mev linear accelerator at the Lewis Research Center. The procedure 
would involve first calculating the energy (velocity) distribution of primary electrons 
that are emitted from the surface of the gamma irradiated material and then synthesiz- 
ing this distribution by using a linear accelerator beam. From electron beam synthesis 
of nuclear reactor irradiation experiments, the extent and method by which energetic 
primary electrons contribute to increased current flow in irradiated diodes could be 
determined. 

The transfer of energy between photons and matter is important in the calculation of 
energetic electron emission from the surface of an irradiated target. In this energy 
transfer process, photon energy is converted to kinetic energy of electrons and to rest 
energy of electron-positron pairs, and then the kinetic energy of the primary electrons 
is dissipated along their tracks because of interactions with atomic electrons and nuclei. 
Electrons produced near the surface, however, may not deposit all of this energy and 
thus be able to escape depending on the spatial location and initial energy of the electron. 
This report is limited to calculating the current density and energy spectra of these es- 
caping primary electrons. 

Gamma ray interactions with matter have been extensively studied, and specifica- 
tions of the various interaction cross sections have been reported. Descriptions of these 
processes and data on energy transfer coefficients are found in works by Fano, Spencer, 
and Berger (ref. 3); Grodstein (ref. 4); Berger (ref. 5); and Evans (ref. 6). 

Electron transport through matter and the basic electron interactions which occur 
during flight are described by Rossi (ref. 7), Bethe and Ashkin (ref. 8), Fano (ref. 9), 
and Birkhoff (ref. 10). Different calculational techniques and theoretical models have 
been applied in order to obtain better agreement between calculated results and experi- 
mental data. Monte Carlo methods by Berger and Seltzer (ref. 11), the phase-space 
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time evolution method by Cordaro and Zucker (ref. 12), and the method of discrete or- 
dinates using the code ANISN by Bartine, Alsmiller, Mynatt, Engle, and Barish (ref. 13) 
represent major efforts in the study of electron transport utilizing high-speed large 
computers. The main contribution of the work presented in this report lies in the for- 
mulation of a simple computational method of estimating electron penetration by using 
available experimental data. 

The numerical results presented are for a tungsten target 0. 33 centimeter thick and 
2. 5 centimeters in diameter irradiated by the fission gamma flux in the core of the NASA 
Lewis Research Center Plum Brook nuclear test reactor. The spatial distribution of 
gamma radiation is prescribed according to specifications of geometry and fuel loading 
of the core. The fission gamma energy spectrum is divided into 13 discrete energy 
groups, each one characterized by a monoenergetic gamma flux having the mean energy 
of the group. The generation of primary electrons is deduced from appropriate energy 
transfer coefficients as a function of gamma energy. The spectrum of primary electron 
energies resulting from gamma interactions is approximated by an average energy within 
each specified gamma group. Empirical electron energy-range relations are used in 
conjunction with the probability of escape and the average path length of escaping elec- 
trons to compute the emission of energetic electrons. The derivation of the latter two 
quantities follow along the lines first suggested by Butler. This approach differs from 
an earlier calculation by Shiomi and Higashimura (ref. 14), who adapted the penetration 
theory of Spencer (ref. 15) and the segment theory of Higashimura (ref. 16) in describ- 
ing electron emission from a plane semi-infinite target results from the passage of a 
normal incidence gamma flux. The energy spectra and flux of emerging electrons from 
a uniform plane slab target of finite thickness are calculated by assuming an isotropic, 
monoenergetic beam of photons for each group. The total emission current density and 
power per unit area are then computed by summing over the 13 discrete energy groups. 


SYMBOLS 

B g binding energy, eV 

c speed of light, m/sec 

-IQ 

e charge of electron, 1.6x10 C 

o 

secondary electron energy flux within target, W/cm 

2 

gamma energy flux, W/cm 

2 

J egc current density of escaping secondary electrons, p A/cm 
K linear pair production attenuation coefficient, cm"'*' 
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m 0 electron rest mass, kg 

P probability of escape for secondary electron 

cbL 

o 

R electron range, g/cm 

r electron path length defined in fig. 1 

2 

T target thickness, g/cm 

2 

W energy flux of escaping secondary electrons, W/cm 

x dimension defined in fig. 1 

initial energy of electron, eV 
e ocn energy of electron upon escape from target, eV 
energy lost by electron within target, eV 
gamma energy, eV 
dimension defined in fig. 1 
critical angle for escape 

linear coefficient of gamma energy absorption, cm - * 

total linear gamma energy attenuation coefficient, cm 

2 

density of target, g/cm 
Compton absorption coefficient, cm - * 

Compton scattering coefficient, cm - * 
photoelectric linear absorption coefficient, cm * 
azimuthal angle 



PROBABILITY OF ESCAPE 

The probability that an electron, generated by photon interaction within a target, 
will escape to the surface is derived in this section. The target is assumed to be an 
infinite slab which is immersed in an isotropic monoenergetic gamma flux field. The 
electrons produced by photoelectric, Compton, or pair production interactions will like- 
wise have an isotropic distribution; however, the electron energy distribution associated 
with a given gamma energy will be simplified by assuming each electron has an energy 
equal to the average. The methods and implications of this averaging process are dis- 
cussed in the section CALCULATION OF EMITTED ELECTRON CURRENT. 
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Since the cross section is large, electron-electron scattering is the major mecha- 
nism of energy degradation of the primary electron as it traverses the medium. It is 
because of the almost continuous loss of energy in multiple electron-electron interactions 
that primary electrons have an almost precise range R. For a target whose thickness 
is greater than the range of the scattered electrons, the escape probability is different 
from that for a target of smaller thickness. The derivation is best begun by reference 
to figure 1. 

The line r in figure 1 does not represent as escape path for an electron generated 
by x\ In traversing the medium to the surface, electrons will lose energy in numerous 
collisions, and as a result the escape path will deviate considerably from a straight line. 

The differential probability of escape dP for an electron scattered between 

C Dv 

x = x' + dx and x = x' in an infinite slab is the fraction of the solid angle of the sphere 
of radius R contained within the cone 9 < 9 c : 


dP 


esc 


4 v 



sin 9 dc p d0 


( 1 ) 


The total probability of escape from one side of the target where the thickness is 
greater than the range is obtained by averaging dP throughout the target: 

6 SC 


p 

esc 

T>R 



( 2 ) 



Figure 1. - Diagram for computing 
probability of escape for electron 
scattered at x 1 in target of thick- 
ness T. 
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When the thickness of the target is less than the range, the probability of escape 
from one of the surfaces may be written 


P 

esc 
T< R "0 


/•T /•£? /*2 tt 

/ / / sin 9 

h ~r\ Jr\ 4ffT 


0 "0 


sin 9 d<p dd dx = - ( 1 - —) 

2 V 2R 


( 3 ) 


AVERAGE PATH LENGTH FOR ESCAPING ELECTRONS 

When the thickness of the target is greater than the electron range R, as shown in 
figure 1, the path length R of electrons escaping from one side of the slab from x, 
where T - x is less than the range, will be a function of the scattering angle 9: 


r = 


T - x 
cos 9 


(4) 


The mean path length of escaping electrons is computed by averaging r over the 
pertinent ranges of x, cp, and 9: 


f T f C 

J np t> r'l *rv 


r = 


T-R 0 


— — - R^ sin 9 dcp d 9 dx 
cos 0 


/ R f 9c f 2 * 

III sin 0 d cp d0 dx 

•'T-R *'0 •'0 


(5) 


where 


0 C = cos 


-1 T - x 


After the appropriate substitutions, the integral (5) can be solved for T > R to ob- 
tain the simple result 


r 


R 

2 


( 6 ) 
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For thin targets (R > T) the integral (5) is equally valid with the modification that 
the limits of integration on x are from 0 to T. 

The average path length of electrons escaping from one side of the target is 


r 


ln* + A 


T 

2 


1-X 

2R 


( 7 ) 


CALCULATION OF EMITTED ELECTRON CURRENT 

The quantities which must be known for the computation of high-energy electron 
emission from a target in a gamma flux field are the gamma energy e , the gamma 
energy flux J , the target density p, the total linear gamma energy attenuation coeffi- 
cient p^., and the linear coefficient of absorption of gamma energy electrons p g . The 
development will be restricted to targets thicker than the average electron range. 

Even a monoenergetic gamma flux field will produce electrons with a wide spectrum 
of energies. In this calculation, the distribution of electron energies is approximated by 
an electron flux of constant energy equal to the average. This first step in the calcula- 
tion is to find the average energy of the electrons resulting from the attenuation of 
the gamma flux in the target. When p g and p t are known, the equation for e g may 
be written 


e 


e 



(8) 


When tabulated values of p g and p^. are not available or when a theoretical analysis of 
the energy absorption process is desired, equation (8) may be expressed (ref. 6) 


€ e = € y 


cr + cr + r + K 
a s 


(9) 


where a is the Compton absorption coefficient, r is the photoelectric linear absorp- 

a, 

tion coefficient, B g the binding energy of the target atoms, K the linear pair production 
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attenuation coefficient, m Q the rest mass of the electron, e the electronic charge, 
c the speed of light, and cr the Compton scattering coefficient. 

The extrapolated range R of the electrons produced as a result of photon inter- 
actions in the target is found in the literature for the appropriate material and energy 
(refs. 17 to 19). The range-energy relation for tungsten taken from reference 17 is 
plotted in figure 2. 



The gamma energy flux drops off exponentially as it penetrates the target material 
(ref. 6). The primary electrons produced do not have as great a penetration distance as 
the gamma rays, so that 


* « — ( 10 ) 
p p e 


Therefore, the power per unit area absorbed by the electrons within a distance 
R/p from the surface may be expressed approximately as 


> e = -'y 


^Ru 


s\ P 


( 11 ) 
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For the thick targets treated in this report the average path length r of escaping 
electrons is shown in equation (6) to be R/2. The average energy loss e L for elec- 
trons can then be found from figure 2. The average energy e of escaping electrons 

G bb 

is then computed from 


€ 


esc 


= e. 


e L 


( 12 ) 


The energy flux W of escaping energetic electrons can now be written 


W = 


Vesc P 


esc 


(13) 


Ideally, the escape of electrons from a depth greater than one range is not possible; 
hence, only that portion of the target which is one range thick need be considered. From 
equation (2) the probability of escape is simply 1/4, so that equation (13) may be written 


W = 


J e 

e esc 

4? 


(14) 


Finally, the escape current density J is simply expressed as 

G bv 


J 


esc 



(15) 


RESULTS FOR TUNGSTEN 

The calculation presented in this section refers to a tungsten target 0. 33 centimeter 
thick placed within the gamma flux field approximately 4 centimeters from the core of 
the NASA Plum Brook Reactor. This system is of interest because of its possible appli- 
cation to a thermionic converter designed to operate within the core of a nuclear reactor. 

The gamma flux distribution at the point of interest was obtained by using the com- 
puter code QAD P5A (refs. 20 and 21) and the gamma energy spectrum of Groshev 
(ref. 22). The gamma energy spectrum of a nuclear reactor is nearly continuous and 
not a monoenergetic beam as called for in the computation; however, by splitting the 
spectrum into 13 energy groups, the average energy and energy flux of each may be 
specified. This is done in table I, second and third columns. The linear attenuation and 
absorption coefficients are obtained from published data (ref. 23) and for convenience 


9 



o 

» 2 - 

c ~ c 

o x <u 

' — a 


co^ioiocD'^'^cocOTfcooC'- 
CM-^^CDCOCMCOt— ‘C'-’^CMO’-* 
1-4 o o o O H H N CM CO ^ lO C*” 




c 

° OJ >, 

>h 2 * bn 

u o £ 

<D 03 c 


o 

fcD -Q Sh 

|ll 


•5 Si j jg 

% s |w s 


S' - c5 .. £ 

m 3 bO v tl) O 

C y b ^ \ 

<U g 5 £ 





CO 

CO 


i— < 

LO 


CO 

co 

05 

LO 

CM 

05 



LO 

1—4 

CO 

LO 

05 

LO 

CO 

05 

CO 

CO 

O 



1 — I 

CM 

CM 

CM 

CM 

CO 

co 

CM 

CM 

1-4 


o 

o 

o 

O 

O 

O 

O 

o 

o 

O 

O 

o 

o 

d 














o 

o 

CM 

t- 

CM 


LO 

CM 

05 

1*4 

LO 

t- 



c— 

O 

1-4 

05 

CO 

CM 



c- 

t- 


o 

1-4 

CM 

o< 

to 

CO 

CO 

o 

CM 

"3< 

CO 

CO 

CO 

d 







1—1 


1-4 



CM 



LO 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 


o 


1 -4 


o 

CO 

CM 

o 

c~ 

LO 

LO 

O 

i-4 

CM 

CO 

LO 

c- 

05 

o 

CO 

CO 

CO 

1-4 


1-4 

d 







1—4 

1—4 


CM 

CM 

CO 

CO 

c— 


rr 

LO 

LO 

LO 

CO 

LO 


to 

CO 

o 


t-4 


05 

CM 

CO 

CO 

CM 

o 

t- 

1-4 

CO 

CO 

o 

rH 

i-4 

r-4 

CM 

CM 

CM 

CO 

co 

CM 

CM 

^4 




ooooLOxrioioiOLOioo 

mC-*-4C0£~-I^C005Ot-<CMC0CD 

OinWWNNNNCOC^COCON 

cooooooooooooo 


C 


















O 


CM 





CM 

o 

LO 

O 

O 

o 

o 

o 

lO 

O 


CD 


£ 





i—i 

1—1 

1—4 

CM 

CO 

CO 

o 

05 

r- 


o 

bl) 

c 

os' 

o 


O 

o 


CM 

CO 


LO 

CO 

CO 

o 


CO 

CO 

3 

rt 

u 


bB 


d 










1—1 

1-4 

1-4 

c 


















o 

>> 






in 

CM 

c- 

CM 

o 

LO 

CM 

05 

1-4 

LO 

I> 


be 


> 






LO 

05 




1-4 

CM 

CM 


"o 

u 

CD 

o 


CM 

CO 

co 

05 

CM 

VO 

05 

CO 

00 

CO 

CO 

CO 


<D 

c 


s 


o 

• 



_4* 

1-1 


CM 

CM 

CO 

CO 


LO 

w 

CD 

















ji 

O 

C 
















c3 

4-4 


c 







LO 

LO 

LO 


LO 

C— 

c- 


fi i 

«—• 

■w 

CD 

=L 



LO 

CO 

CO 

CO 

05 

CM 

CO 

CO 

o 

CM 

CO 

O 

O 

o 

CD 

o 

c 

o 

\ 

<D 

3. 

CO 

d 

LO 

LO 

LO 

CO 

CO 

CO 

c— 


c— 

CO 

CO 

00 

■4 

^4 


CD 















£ 

o 

— ’ 

O 















Ph 

tfi 


o 
















S W) - P 

s S 3 

rt c r: ' K x 

O u £ 


(MONHCO^iMtNCO^^O© 

00^^-00505^0100510^0 

LOi-4i^tr-ooioo5coocoiocot- 

HO)^N0O^to'cNJCOifiCOiri 


r- >> ^ 

c bJ) „ > 

£ S-. - QJ 

5 <s s 

O <1) 



hcociN^^win co 

COHHt-OOOCn^lN 

c\j q£) •—* lo 05 co tr~* cm t' - cm 

o i-4 »— < CM CM CO CO 


^CMCo^mcot-cocoo^ooco 


10 


6. 334 




























Figure 4. - Absorpt 



0 



plotted in figures 3 and 4. The calculation then proceeds as described in the preceding 
section. 

The results of the calculation, the power per unit area and current of escaping elec- 
trons, are shown in figures 5 and 6. The totals for all energy ranges are 0. 3109 watt 
per square centimeter and 0.4506 microampere per square centimeter. By comparison, 
in a cesium -filled thermionic converter operating in the power producing mode, the loss 
in the plasma is of the order of 5 watt per square centimeter. 

This value for the current density of electrons in the million volt energy range has 
been utilized in experiments in which the effect of gamma radiation is simulated by the 
electron beam of a linear accelerator. The results of some of these experiments were 
reported by Bacigalupi (ref. 24) and Forman (ref. 2). In tests of a xenon -filled thermi- 
onic diode, excellent agreement was obtained for current, voltage curves taken in a nu- 
clear reactor and in a linear accelerator which was delivering a current density pre- 
dicted by this method of calculation. 

In addition to the primary electrons having energies of the order of 1 MeV, second- 
ary electrons of energy less than 100 volts will also be emitted. These secondary elec- 
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Normalized energy flux of escaping electrons, 
(WA/J x 100, percent 
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tr.ons have been found to be more numerous than the primary electrons by some inves- 
tigators (refs. 25 and 26). The effect of the low-energy secondary electrons presumably 
could in some instances be as important as that of the primary electrons in altering 
thermionic converter performance. 

Of more general interest may be the plot of W /j shown in figure 7. Within the 
accuracy of the computation, figure 7 should apply to any tungsten target irradiated by 
a gamma flux field of any specified energy spectrum as long as the thickness of the tar- 
get is greater than or equal to the range. 

The rise in W/.X for incident gamma rays above 0. 7 MeV can be readily explained. 
If equations (11) and (12) are substituted into equation (14), an expression is obtained for 

W//: 



As the gamma energy is increased, the range increases and the ratio of electron 
energy lost to original electron energy decreases, both factors tending to increase 
W /j^. The absorption coefficient, figure 4, drops rapidly up to about 2 MeV and then 
increases slowly for higher values of and thus accounts for the shape of figure 7. 

It should be emphasized once more that figure 7 applies only to the case where 
range does not exceed target thickness. Therefore, for any target of finite thickness 
the ratio of W to / will not increase indefinitely for incident gammas of higher 
energy. 


SUMMARY OF RESULTS 

This report presents a method for estimating the emission of high-energy electrons 
from a metal target irradiated in a gamma flux field. The calculation was based on 
photon energy-transfer coefficients, the electron energy range relations, the probability 
of electron escape, and the average path length of emitted electrons. The latter two 
quantities were derived in this report. 

The magnitude of electron current density and energy flux resulting from the place- 
ment of a tungsten target in the gamma radiation field of a nuclear reactor was com- 
puted. In this numerical example, the continuous gamma energy spectrum was broken 
into 13 monoenergetic groups, and the photon, electron interactions were treated for 
each group separately. For a 0. 33-centimeter-thick tungsten target located 4 centime- 
ters from the core of the NASA Plum Brook Reactor, a total primary current density of 
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0.4506 microampere per square centimeter and an energy flux of 0. 3109 watt per square 
centimeter were estimated. 

This calculation has been used successfully to predict a starting point for the simu- 
lation with the beam of a linear accelerator an effect observed in a nuclear reactor. The 
agreement between reactor and linear accelerator data reported so far may be fortuitous, 
considering the many approximations on which this calculation was based. The results 
obtained should be taken only as an estimate of the million-volt -electron emission caused 
by the gamma radiation. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 26, 1972, 

503-25. 
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"TAe aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 
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